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Abstract. Two super-integrable and super-separable classical systems which can be con- 
sidered as deformations of the harmonic oscillator and the Smorodinsky-Winternitz in two 
dimensions are studied and identified with motions in spaces of constant curvature, the 
deformation parameter being related with the curvature. In this sense these systems are 
' to be considered as a harmonic oscillator and a Smorodinsky-Winternitz system in such 

, bi-dimensional spaces of constant curvature. The quantization of the first system will be 

carried out and it is shown that it is super-solvable in the sense that the Schrodinger equa- 
, tion reduces, in three different coordinate systems, to two separate equations involving only 

' one degree of freedom. 

l> ! 
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1 Super-integrable systems 

\^ • There are few integrable systems in the Arnold-Liouville sense: Hamiltonian systems in a 2n- 

. dimensional symplectic manifold for which there exist n functionally independent constants of 

5^ \ motion /j in involution (including the Hamiltonian H itself), i.e. 

{/^,/i}=0, Vi,j = l,...,n, d/iA---Ad/„^0. 

The system is said to be super-integrable when it is integrable and there exists a set oi m > n 
functionally independent constants of motion, i.e. it possesses more independent first integrals 
than degrees of freedom. The existence of these additional first integrals gives rise to a higher 
degree of regularity in the phase space (e.g. there exist periodic orbits) since the trajectories are 
restricted to submanifolds of dimension lower than n. In particular, a system with n degrees of 
freedom possessing 2n — 1 independent first integrals is said to be maximally super-integrable. 

It is also well-known [1] that the only central potentials in which all bounded orbits are 
closed (periodic) are given hy V = (l/2)u;Qr^ and V = —k/r. From a modern point of view the 
existence of closed trajectories is considered as a consequence of the existence of the maximal 
number of functionally independent integrals of motion; thus, the result obtained by Bertrand 
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is a proof of the super-integrability of these two systems. In fact, the harmonic oscillator 
shares with the Kepler problem a very distinguished property not only in Classical but also 
in Quantum Mechanics, where super-integrability leads to energy levels depending on a single 
quantum number. 

Fris et al. studied in 1965 the Euclidean n = 2 systems which admit separability in two 
different coordinate systems [2], and obtained four families 14, r = a, 6, c, d, of super- integrable 
potentials with constants of motion linear or quadratic in the velocities (momenta). The two 
first families 

ya = ^ujQ{x +y ) + — + —, Vb = -ujQ{Ax +y ) + k2X + 
z X y I y 

can be considered as the more general Euclidean deformations (with strengths k2, k^) of the 1 : 1 
and 2 : 1 harmonic oscillators preserving quadratic super-integrability (the other two families, Vc 
and Vd, were related with the Kepler problem). The super-integrability of Va, which is known as 
the 'Smorodinsky-Winternitz' (S-W) potential, has been studied by Evans [31 0] for the general 
case of n degrees of freedom. 

A large number of papers have been published on super-integrability in these last years, 
most of them related with quadratic superintegrabilty (see [5] for the proceedings of a workshop 
on super-integrability and Refs. [HI [3 El [9] for some very recent studies on super-integrable 
systems with integrals quadratic in momenta). The idea is that if we call super-separable 
a system that admits Hamilton-Jacobi separation of variables (Schrodinger in the quantum case) 
in more than one coordinate system, then quadratic super-integrability (i.e., super-integrability 
with linear or quadratic constants of motion) can be considered as a property arising from 
super-separability. We note that these studies also include non-Euclidean Hamiltonian systems 
[lilllliniliaiTllllIlllSinillBllISEOlE^ and that in both cases, Euclidean and non- 
Euclidean, many of these systems are closely related with the harmonic oscillator. 

The rareness of integrable systems leads to the following question. Is it possible to deform 
a super-integrable system but preserve super-integrability of the system? We report here some 
results of previous works where such question is affirmatively answered. We shall show that this 
is possible for two important examples: the isotropic bi-dimensional harmonic oscillator and the 
Smorodinsky-Winternitz system. 

A very interesting example of a one-dimensional nonlinear oscillator depending of a parameter 
A was studied in 1974 by Mathews and Lakshmanan [231 [24] and it has recently been proved [25] 
that this particular nonlinear system can be generalized to the two-dimensional case, and even 
to the n-dimensional case: these higher dimensional systems are super-integrable deformations 
admitting 2n — 1 quadratic constants of motion. It is also proved that there is a related A-de- 
pendent two-dimensional oscillator isotonic oscillator that is a super-integrable deformation, i.e. 
a A-dependent version of the Smorodinski-Winternitz system [261127] . Actually, the deformation 
introduced by the parameter A modifies the Hamilton-Jacobi equation but preserves the exis- 
tence of a multiple separability. Moreover, we point out that a geometric interpretation of the 
higher-dimensional systems was proposed in relation with the dynamics on spaces of constant 
curvature. 

This paper must be considered as a survey summarizing recent works by the authors on 
properties related with the integrability and super-integrability of certain two-dimensional A-de- 
pendent systems related with the harmonic oscillator and it is neither a review of the whole field 
nor a comparison with other approaches of different authors. It is mainly focused on the study of 
deformations of super-integrable systems that do not alter the super-integrability structure. It 
starts with a brief review of the basic properties of the above mentioned isotonic oscillator, S-W 
system and Mathews-Lakshmanan oscillator and some two-dimensional A-dependent oscillator- 
like systems, first in in the classical approach and afterwards in their quantum counterparts. 
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Roughly speaking this survey is divided into three main parts: first, the analysis of the 
classical A-dependent oscillator including the existence of Hamilton-Jacobi multiple separability 
(existence of alternative coordinate systems in which the corresponding A-dependent Hamilton- 
Jacobi equation separates). Second, the relation of this A-dependent nonlinear model with the 
harmonic oscillator on the three spaces of constant curvature (5^,E^,i?^). The final part is 
devoted to the analysis of these systems from the quantum viewpoint. 

The second part presents a geometrical approach and proves that these systems can be 
considered in two different ways: either as a nonlinear deformation of a linear system, or simply 
as a model of the oscillator on spaces of constant curvature. In the first case the parameter A 
represents the strength of the deformation and in the second one the curvature of the space. 
Finally, in the third part, devoted to the quantum version of this A-dependent oscillator, it 
is proved that the system is exactly solvable: the Hamiltonian can be factorized and the wave 
functions and energies can be explicitly obtained. It includes some points such as: (i) Analysis of 
the transition from the classical A-dependent system to the quantum one, (ii) Exact resolution 
of the A-dependent Schrodinger equation, factorization method and existence of operators A 
and A~^, and shape- invariance property, (iii) Schrodinger multiple separability and quantum 
super-integrability. 

In more detail the structure of the paper is as follows: In next section we describe shortly two 
one-dimensional classical systems which can be considered as generalizations of the harmonic 
oscillator and the isotonic oscillator: the A-deformed nonlinear oscillator and the deformed 
isotonic oscillator. Section 3 deals with two-dimensional generalizations of these classical systems 
and their deformations and Section 4 studies the separability properties of their Hamilton Jacobi 
equations. A geometric interpretation is given in Section 5 and a quantization of the nonlinear 
oscillator, is carried out in Section 6, and the quantum spectrum s computed by using the 
traditional power series expansion. The system is shown in Section 7 to admit a shape invariant 
factorization which allows us to explicitly compute the spectrum in an algebraic alternative way. 
Finally, in Section 8 we sketch the method to be used for the corresponding two-dimensional 
system and prove that he Hamiltonian can be written as a sum of three different terms such 
that each one commutes with the sum of the other two, what provides us alternative complete 
sets of compatible observables. 



2 Some one-dimensional classical systems 

We start this section by reviewing the basic properties of two simple classical systems. 



2.1 The harmonic oscillator 

The dynamics of the classical harmonic oscillator in one dimension is given by 

dx dv 2 

— = V, — = —UJ X 
dt dt 

and is described by a Lagrangian L = (1/2) (?;^ — w^x^). A complex variable z = lox + iv may 
be introduced and the equations of the motion become 

dz 

whose general solution \s z = zoe"*'^*, i.e. 

X = xo cos ujt sinwt = Acos{ujt + (p), 

LO 

and therefore the solutions are periodic with angular frequency w, while A and ip are arbitrary. 
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2.2 The isotonic harmonic oscillator 

The isotonic oscillator is described by the Lagrangian [281 ES] 

L = -t; - -ax fc > 0, 

2^2 x^ 

(like a harmonic oscillator plus a centripetal barrier). This is an important example of an 
isochronous system and actually they are the only two rational potentials giving rise to isochro- 
nous systems [30l [31] . 

The Euler-Lagrange equation 

x + a X -\ — ^ = 0, c = —2k, 
x-^ 

is a particular case of the so-called Pinney-Ermakov equation [321 [33] whose general solution 
can be written as x = {\/{aA)\/ {a^A^ + c) sin^(at + (f) — c. 

The corresponding quantum system admits a shape-invariant factorization [M] and therefore 
is solvable by means of algebraic methods. 



2.3 A 1-dimensional nonlinear oscillator 

In 1974 Mathews and Lakshmanan [23] studied the equation of motion 

{l + \x^)x-\xx^ + a^x = Q, A>0. (1) 

The general solution takes the form x = Asm.{u)t + (j)), with the following additional restriction 
linking angular frequency u and amplitude A: 

2 

UJ — 



1 + ■ 

The equation ([T|) is the Euler-Lagrange equation for the Lagrangian: 
1 1 

It describes a system with nonlinear oscillations with an amplitude dependent frequency (or 
period) . We can also allow negative values for A [2^ , but when A < the values of x are limited 
by the condition \x\ < l/y^|A[. In the limit A — > we recover the equation of motion and 
both the Lagrangian of the harmonic oscillator and the frequency become independent of the 
amplitude. Therefore, the system can be seen as a deformation of the harmonic oscillator. 

It can also be seen as an oscillator with a position-dependent effective mass which depends 
on A: 

1 

mx 



1 + Xx'^' 

The Hamiltonian for such a system is: 

i^A(x,p) = i(l + A.V + ^Y^- 

The important fact is that there is an interesting generalization to n = 2 or even arbitrary n [25j . 
Note that the Lagrangian for the one-dimensional free-particle (i.e. for a = 0) 

i.(x.»J = T,(A) = i^ 
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is invariant under the vector field 



Xi{\) = v/l + Ax2— + 



d Xxvr d 



which is the natural lift to the phase space M x M of the vector field in M 

X,{\) = v/l + Ax2 — , (2) 
i.e. X*(A)(ri(A))=0. 
2.4 A deformed isotonic oscillator 

We study next a deformed isotonic oscillator in one dimension described by the Lagrangian [26] : 
1 1 a^x"^ k 

Lx{x,v„k) - ^Y^x^ - 2TTX^ ~ ^' 

Here A can be any real number, but when A < the possible values of x are such that |x| < 
The Euler-Lagrange equation is 

(Px Xx f dx\^ o?x 1 + Ax^ 

dt^ ~ 1 + Ax2 V j ^ 1 + Ax2 " x3 " ' 

and one can see that the general solution for bounded motions is |26j : 

X = Y^(w2A4 - 2k) sin2(u;t + (/>) + 2k, 
where 

ujRi [2k + (w^yl^ - 2k) sin2(a;t + ,/-)]^ = 

with 

Ri = Xlv^A'^ - (a2 - c^2 _ 2A;A2)a2 + 2A:A. 

It can be done similarly for unbounded motions, the functions s'm{Lvt + (j)) being then replaced 
by functions sinh(r2t + <p). There also exist limit unbounded motions of the form 



x = ^/{AtTW+C. 

It will be shown that as in the harmonic oscillator case, there is an interesting generalization 
to the n = 2 case. 

3 Some two-dimensional classical systems 
3.1 The harmonic oscillator 

The Hamiltonian of a two-dimensional classical harmonic oscillator is (for simplicity m = 1): 
H{x,y,pcc,Py) = ^{pI +pI) + ^{ujfx'^ + ujjy'^) , 
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and we can easily check that the energy functions for each degree of freedom are constants of 
motion: 

Ii = E^ = ^ {pI + ujjx^) , I^ = Ey = ^{pl + uly^) . 

The rational case, for which wi = riiujQ, uj2 = n2UJo, with ni,n2 £ N, is super-integrable j35]. In 
fact, let Kx and Ky be defined by = Px + iniuox, Ky = Py + in2UJoy- The Hamiltonian H 
and the canonical symplectic form can be expressed in terms of such functions as follows: 

H = l [KxK + KyK) , f^o = TT^dKx A dKl + -^dKy A dK* 

and therefore, as the fundamental Poisson brackets are {Kx,K*} = 2iniujQ and {Ky,Ky} = 
2in2UJo, the evolution equations are 

= iniLOoKx, = -in2LOoK*. 

Hence, the complex function J defined as Jni,n2 = ^x^i^y)^^ is a (complex) constant of 
motion which determines two different real first integrals: Js = Im(J„j^„2) and I4, = Re(J„^^„2)- 
They are polynomials in the momenta of degrees ni + n2 — 1 and ni + 712, respectively. 

Two particularly important examples are the isotropic harmonic oscillator (for ni = n2 = 1) 
and the case ni = 1, 77-2 = 2. In these cases the four constants are not independent but only 
three of them are independent. For instance, in the isotropic case the constants of motion take 
the form = PxPy + uj'^xy and I4, = xpy — ypx- 

3.2 The Smorodinsky— Winternitz system 

A 2-dimensional generalization of the isotonic oscillator with rotational symmetry would be 



Vr^ = -UjI{x'^ + y"^) + 



There is however another super-integrable generalization given by the Smorodinsky- Winternitz 
potential [2] 

Vsw = 2^o(a^ +y ) + ^ + ^' 

which is a set of two non-interacting isotonic systems with the same frequency and in general 
different constants k2 7^ A;3. Of course, the case A;2 = ^3 = reduces to the usual isotropic 
harmonic oscillator. The energies of each degree of freedom are constants of the motion and 
there exists a third constant of motion given by 
2 ^2 

C = J"^ + k2— - k^ — -, with J = xvy - yvx- 
X y 

3.3 A nonlinear oscillator in two dimensions 

It has recently been proved [25] that there exists a 1-parameter dependent generalization of the 
nonlinear oscillator for the 2-dimensional case with the following requirements: 

1. The kinetic term T2(A) is a quadratic function of the velocities invariant under rotations. 

2. T2(A) is invariant under (the tangent lifts of) the vector fields ^i(A) and X2{X) given by 

Xi(A) = v/TTa^I-, X2{X) = Vl + Ar2|-, 
ox ay 

which are extensions to of the vector field Xx{X) in the n = 1 case given by ([2]). 
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These conditions lead to the following form for the kinetic energy: 

1 1 

2 1 + Ar2 



The term XixVy — yvx)"^ represents a two-dimensional contribution that cannot appear in the 
one-dimensional case. When A < this function will have a singularity at 1 — jAjr^ = and we 
should restrict our dynamics to the interior of the circle + y'^ < 1/|A| where T2(A) is positive 
definite. 

The kinetic energy T2(A) is determined by the A-dependent metric 

ds^{X) = + \y^)dx^ + {l + \x^)dy^ - 2\xydxdy]. 

T2(A) remains invariant under the actions of the lifts of the vector fields Xi(A), X2{X), and Xj, 
given by 

Xi(A) = Vl + Ar^A, X2(A) = Vl + Ar2|-, Xj = x^-y^, 
ox oy oy ox 

namely 

OX Vl + Ar2 dv^ dy Vl + Xr^ ovy 

X* ^ ^ 9 d 

^ ^ dy ^ dx dvy dv^ ' 

These vector fields close on a Lie algebra: 

[Xi(A),X2(A)] = XXj, [Xi(A),Xj] =X2(A), [X2{\),Xj] = -Xi(A), 

which is either isomorphic to so(3,M) when A > 0, to 5o(2, 1) when A < 0, or to the Lie algebra 
of the Euclidean group for A = 0. 

The appropriate generalization for the potential of the nonlinear two-dimensional A-dependent 
oscillator is given by 

V.{x,y)-^' + 



2 l + A(x2 + y2) 



This bi-dimensional nonlinear oscillator is completely integrable [25j . because one can show 
that, if Ki and K2 are the functions 

Ki = Pi (A) + ia , ■ K2 = P2{\) +ia ^ 



Vl + Ar2' ' Vl + Ar2' 

with 

-^i^^) = ^r=^' ^2(A) = ^====, J = xvy-yv^, 
V 1 + Ar^ V 1 + Ar^ 

then the complex functions Kij defined as Kij = KiKj, i,j = 1, 2, are constants of motion. 
In fact, the time-evolution of the functions Ki and K2 is 

dt^' - TTxT^^'' dt^' - TTx^^^'' 

from which we see that the complex functions Kij are constants of the motion. 

Therefore the system is super-integrable with the following first integrals of motion: 

/i(A) = |J^i|2, /2(A) = |iC2|', h = Im{Ku) = a{xVy-yv^). 
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3.4 Nonlinear Smorodinsky— Winternitz system 

We can consider a deformed Smorodinsky-Winternitz system with Lagrangian [26] 

T ^ r 2 , 2 , ^/ ^21 x2 +y2 k2 h 

2(1 + Ar^j '- -'2 1 + A(2;^ + y^j 

The E-L equations cannot be directly solved in a simple way, but one can check by direct 
computation the existence of three independent constants of motion: 

1 + A(x^ + y^j x"^ 

This proves that for any value of A this system is completely integrable. Note that for 
A = the previous constants of motion reduce to the three first integrals of the S-W potential. 
Consequently, this system is a deformation of the S-W potential which preserves the completely 
integrability. 

4 Hamilton— Jacobi super-separability of the previous systems 

The Legendre transformation for a Lagrangian Lx = T2{X) — V{x) is given by 
(1 + \y'^)vx - Xxyvy _ (1 + ^x^)vy - Xxyv^ 



~ 1 , \^2 ' Py 



1 + Ar2 ' 1 + Ar2 

Note that xpy — ypx = xvy — yvx ■ 

The general expression for a A-dependent Hamiltonian is 

H{X) = ^[pl+pl + X{xpx + yPyf] +^aV(x,y), 
and hence the associated Hamilton-Jacobi equation is 
dSy , /dsV . . f dS , dsV 



ox J \oy J \ ox ay J 

This equation is not separable in (x, y) coordinates but there exist three particular orthogo- 
nal coordinate systems, and three particular families of associated potentials, for which such 
a Hamiltonian admits Hamilton-Jacobi separability: 

(i) In terms of the new coordinates {zx,y), Zx = x/ \/l + Xy^, the Hamilton-Jacobi equation 
becomes: 

so if the potential V{x, y) can be written in the form 

then the equation becomes separable. 

The potential is therefore integrable with the following two quadratic integrals of motion 

Ii(A) = {l + Xr^)pl + a^Wi{zx), 
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72(A) = (1 + Xr^)pl - AJ2 + (w2{y) - Yj^Wi{z^)^ . 

and Hx can be written as: Hx = (l/2)(/i(A) + /2(A)). 

(ii) Similarly, using coordinates {x,Zy), Zy = y/\/l + Ax^, the Hamilton- J acobi equation 
becomes: 



(1 + Xx^) 



2^2(^S 

dx J \ dz. 



+ 1 + Az, 



dS_ 

V 



^ +a^{l + Xx^)V = 2{1 + Xx'^)E, 



therefore, if the potential V{x, y) can be written on the form 

W2{Zy) 



V = Wi{x) + 



1 + Ax2 



then the equation becomes separable. 

The potential is therefore integrable with the following two quadratic integrals of motion 



Ax^ 



/i(A) = (1 + Xr^)pl - AJ2 + [Wi{x) ^ ^ 
h{X) = {l + Xr^)pl + a^W2{zy) 



Wl{Zy) , 



and Hx can be written as: Hx = (l/2)(/i(A) + /2(A)). 
(iii) In polar coordinates (r, (j)) the Hamiltonian H{X) is 



so that the Hamilton-Jacobi equation is given by 
2 1 / oo\ 2 



V ' \dr J \ d<p 

If the potential V is of the form 



aV(r,0) = 2E. 



then the equation admits separability 

2 



+ r^(a'^F{r) - 2E) + 



(-)' 



+ a'^G{(j)) = 0. 



The potential V is integrable with the following two quadratic integrals of motion: 



Ii(A) = (1 + Ar2)p2 + ^_p2 ^ 



a 



1 - 



F{r) + -^GicP) 



l2{X)=pl + a^G{4>) 
and Hx can be written as: Hx = (l/2)(/i(A) + /2(A)). 
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4.1 Super-separability of the nonlinear oscillator 

Using the previous results one can see that the potential 



a 



2 l + A(x2 + 2/2) 



is actually super-separable (separable in more than one system of coordinates) [25]. In fact, the 
potential for the deformed oscillator we are considering can be alternatively written as follows: 



a 



1 



1 + Az2 



2 1 + Ay2 

Consequently, the Hamiltonian 
1 



a 



1 



2 1 + Ax2 



1 + Az2 



q2 r2 
Y 1 + Ar2 ■ 



admits the following decomposition 

where the three partial functions and are: 



H2{X) = \ 



;i + Ar2)p2+„2 



1 + Ar2 



'l + Xr')pl + a\y 



H-i = ^{xpy - ypxf, 



+ Ar2 

and each one has a vanishing Poisson bracket with 

{H,Hi{\)]=Q, {F,//2(A)} =0, {F,F3(A)}=0 



4.2 Super-separability of the deformed S-W system 

The important point is that, in a similar manner, the A-deformed Smorodinsky-Winternitz 
potential 



_o? ( + j/^ \ k2 h 

2 Vl + A(x2 + y2) y + ^2 + y2 



can be alternatively written in the following three different ways [26] 



a 



X,k 



1 



2 Vl + Ay2 
2 / 1 



a 

Y 

2 



1 + Ax2 

r2 
1 + Ar2 



1 + Az2 



+ 



1 



1 + Ay2 ; 2:2 



+ 



1 + Az2 



+ 4 + 



1 + Ax2 J ^2 



+ 



r2 cos2 I 



r2 sin^ I 



Therefore, it is super-separable since it is separable in three different systems of coordinates, 
{zx,y), {x,Zy), and (r, (/>). This remarkable property means that the Hamiltonian 

Hx,k = \[pI+pI + K^Px + yPyf] + VA,fc 
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admits the following decomposition Hx,k 
tions Hi, H2, and H3, are given by 



Hr, 



H, 



Py 



[1 + Xr^)pl + 



a 



a 



1 + Ar2 

,2 



Hp^ + — XHj, where the three partial func- 
1 + Ay2' 



y 



1 + Ar2 
2 



+ k2 



1 + Xx" 

y2 



1 2 y -^^ 

Hj = {xPy - yPx) +k2^ + /C3-2 • 



Each one of these three terms has a vanishing Poisson bracket with H for any value of the 
parameters A, k2 and fcs 

{-ffA,fc)-f^P;r} = 0, [Hx^kiHpy] = 0, {Hx,k,Hj] = 0. 

Consequently, the Hamiltonian can be written as a sum, not of two, but of three integrals of 
motion. The third one represents the contribution of the angular momentum J to if, with the 
parameter A as a coefficient, therefore vanishing in the limit A — > 0. 



5 A geometric interpretation 

The existence of additional constants of motion for the harmonic oscillator in a spherical ge- 
ometry was studied by Higgs in 1979 [SB]. The Higgs approach considers the motion on S^, 
embedded in the Euclidean space E""*"^, by means of a central (also known as gnomonic) projec- 
tion on a plane IT" tangent to 5"" at a chosen point. In this way he could study the properties 
of the spherical version of the Pradkin tensor. A different alternative approach is discussed 
in [TTl [20] by using curvature-dependent trigonometric and hyperbolic functions, where use is 
made of the curvature k as a parameter so that the dynamics can be studied at the same time 
in the sphere S*^ and in the hyperbolic plane H^. 

In differential geometric terms, the three spaces with constant curvature, sphere S'^, Euclidean 
plane E^, and hyperbolic plane H^, can be considered as three different situations inside a family 
of Riemannian manifolds = {S'^,¥? , H'^) with the curvature ac as a parameter k G M. In 
order to obtain mathematical expressions valid for all the values of k, it is convenient to make 
use of the following K-trigonometric functions 

cos ^/kx if K > 0, 
C^{x) = { 1 if K = 0, Sk(x) = 

COsh-v/ — KX if K < 0, 



and the K-dependent tangent function Tk(x) defined in the natural way, Tk{x) = Sk(x)/ Ck{x). 
The fundamental properties of these curvature-dependent trigonometric functions are 




and 



C'(x) + kS'(x) = 1, 



Ck(2x) = C^(x) -kS^(x), ^Sk(x) = Ck(x), 

S«(2x) =2Sk(x)Ck(x), ^a(x) = -kSk(x). 

ax 
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If in the Lagrangian of the nonlinear oscillator in one dimension 

we consider the change of variable x = Sk{u), where A = — k, we find that 

1 + Ax^ = 1 - = 1 - K Sliu) = diu) 

and as dx/du = Ck{u), we see that 
dx 

X = —u = Ck[U)U 
du 

and therefore the Lagrangian becomes 

the vector field Xx turns out to be 

^ du d d 

dx du du ' 

and the metric giving rise to the Lagrangian g = du (iS> du. We note that the expression T^('u) 
of the above one-dimensional potential appears as directly related with the potential of the 
two-dimensional harmonic oscillator on a space with curvature k. 

The expression of the differential element of distance in geodesic polar coordinates {p, (p) 
on the family = (S^, E^, iJ^), can be written as follows dsf. = dp^ + si,{p)d<l)^ , so that it 
reduces to 

ds\ = dp^ + ( sin^ p)d(jP', dsQ = dp^ + p'^dc/)'^, ds^i = dp^ + ( sinh^ p)d(jP' , 

in the three particular cases n = 1,0,-1 of the unit sphere, Euclidean plane, and 'unit' 
Lobachewski plane, respectively. Note that p denotes the distance along a geodesic on the 
manifold M^; for example, in the spherical k > case, p is the distance of the point to the 
origin (e.g., the North pole) along a maximum circle. 

Therefore, the Lagrangian for the geodesic (free) motion on the spaces (S'^,E^,iJ^) is 

L^{p,(l),Vp,v^) =T,,{p,(f>,Vp,v^) = ^{vl + Sl{p)vl) , 
and the Lagrangian for a general mechanical system (Riemannian metric minus a potential) is 

Lk{p, (t>, Vp,v^) = ^{vf, + Sl{p)vl) - U (p, (j), k). 

The spherical and hyperbolic harmonic oscillators are characterized by the following Lagrangians 
with curvature n [T71 [20] 

L^{p, Vp, v^) = \{vl + Sl{p)vl) - Tl{p), 

i.e. the harmonic oscillators on the unit sphere (Higgs oscillator), on the Euclidean plane, or on 
the unit Lobachewski plane, are: 

Ui{p) = ^ujI tan^ p, Uoip) = U^i{p) = tanh^ p. 



A Super-Integrable Two-Dimensional Non-Linear Oscillator 



13 



Next we study the behavior of under two different changes of variables. 
1. If we consider the ^-dependent change (/o, (^) {t' -,(!>) given by r' = Tk{p), then the 
Lagrangian becomes 



1 / Vr'' r'\l \ 1 



2 /2 

-a r . 



2 Wl + Kr'2)2 (l+Kr'2) / 2 

This function coincides, in the spherical k > case, with the Lagrangian studied by Higgs in |36j . 
In Cartesian coordinates (x, y) it reduces to 

T f \ ^ ^ \ 2 , 2 , I \21-'-2/2 12 2,2 

LHK{x,y,V:,:,Vy) = [v^ + V + K{xVy - yv^) \ - -a r , r =x +y. 

^ (1 + Kr ) ^ 

2. Let us consider the K-dependent change (/o, 4>) i^, 4') given by r = Sk(p), A = —k. Then 
the Lagrangian becomes 

1 / V?: n n\ / 



r / i N 1 2 2A « / 



2 Vl + AH 

Therefore, if we change to Cartesian coordinates {x, y) we arrive to 

1 / 1 \ o? ( 

Lx{x,y,Vx,Vy) = - ( Y^T)^ ) [^x + + A(a;uy - yva;)^] 



2 Vl + Ar2 



2 2,2 

r = X + y . 



This function is just the Lagrangian obtained in [25] as the natural generalization of the one- 
dimensional Lagrangian Lx{x,Vx) for the nonlinear equation ([1]) of Mathews and Lakshmanan. 

We thus have three different and alternative ways of describing the harmonic oscillator on 
spaces of constant curvature: the original K-dependent trigonometric (hyperbolic) Lagrangian 
and the two other approaches, Lhk and Lx, obtained from it. 

The Higgs approach [36] has been studied by many authors (see e.g. [371 EH] and refe- 
rences therein) mainly in relation with the theory of dynamical symmetries. Concerning the 
A-dependent Lagrangian Lx, it has similarities with Lhk but it does not coincide with it. In the 
model of Higgs k (or A) is present in the kinetic term T in a different way and the potential V 
appears as K-independent; this affects to the Hamiltonian formalism. On the other side each one 
of these three formalisms can be used for study of the At-dependent version of the S-W system. 
In the language of the potential, that was studied in [131 [20], is given by 

Ua{r, 4, k) = kiUl + k2Ul + hiUl + fco, ki = \<Jl, 



(Sk(/9)cos(/.)2' (SK(/9)sin, 
with integrals of motion given by 



h{K) = Pi {k)+ujI (t« {p) cos 4>) + 



2ko 



'Tn{p) cos (Pf 



h{K) = P|(k) +a;^(TK(p)sin<A)' + 



(Tk(p) sin( 



T f \ t2/ N , 2/^2 , 2/^3 
= J (/«) + zz::27. + 



cos^ (p sm 
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with -Pi(k), P2{i^) and J{k) given by 

Pi{tv) = {cos (j))vp - {Ck{p) SK{p)sm(j))v^, 

P2{k) = {sin (j))vp + {Ck{p) Snip) cos (j))v^, J{k) = sl{p)v^. 

This S-W system has been studied, using a different approach, in the two- and three- dimensional 
sphere in [10], in the two-dimensional hyperboloid in [12] and in the complex two-sphere in [16] . 

One of the advantages of the Lx approach is that the Euler-Lagrange equations can be directly 
solved and the general solution has a rather simple form that can be interpreted as "quasi- 
harmonic" nonlinear oscillations [SH]; other important advantage is that it is very appropriate 
for the study of the quantum oscillator. 

In what follows we will focus our attention on the quantum Hamiltonian dynamics determined 
by the A-dependent Lagrangian Lx- 

6 The one-dimensional quantum nonlinear oscillator 

Let us consider the quantum case for n = 1. The problem is to define the quantum operator 
defining the Hamiltonian of this position-dependent mass system, because the mass and the 
momentum P do not commute and this fact gives rise to an ambiguity in the ordering of factors. 

Instead of using traditional procedures as Weyl ordering we shall develop an alternative 
method of quantization [39]. We first remark that the vector field 

X,{X) = Vl + Ax2^, 

which was a Killing vector for the metric corresponding to the kinetic energy 

g = {l + Ax^) ^dx (8) dx 

and generates the translations in this Riemann space, does not leave invariant the natural 
measure in the real line but the only invariant measures are the multiples of 

dp= {1 + Xx^y^^^dx. 

This suggests us to consider the Hilbert space d^) and the remarkable fact is that the 

adjoint of the differential operator Vl + Xx'^d/dx in such space is the opposite of such operator. 
The Legendre transformation corresponding to this deformed kinetic energy is defined by 



1 + Ax2 

and then the Hamiltonian function for the free particle is given by 

H ={l + Ax^)^ = i (^^/i + Xx^pY . 

The usual prescription of canonical quantization does not present any ambiguity because the 
linear operator (we put h = 1) 

P = -iVl + Xx^— 
ox 

is self-adjoint in the space £2(M, dfi). Note that with the above mentioned change of coordinates, 
X = Sk{u), dp = du and 

is self-adjoint in the space /^^(M, du), and then P generates translations in these coordinates u. 
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The quantum Hamiltonian operator of the free particle is 

2 



2 \ d 



and in presence of an interaction Vi (x) the Hamiltonian will be 
1 , 1 . d 



We are interested in the case of the nonlinear oscillator for which 
Hi 



1 rP 



1 , d 1 a^x^ 
:Xx— + 



2 dx 2 1 + Ax 



2 ■ 



(3) 



(4) 



If m and h are taken into account, and with the change of parameter, to be more clear later, 
given by 

\ m 

it is enough to define the dimensionless variables 

h 



mB 

and then Hi is written as 
Hi = I3h 



■m(3 



X 



1 /, , '?\ d"^ 1 . 1 .X 



dy^ 2 dy 
The time-independent Schrodinger equation 

Hi<i/ = E = h(58 

becomes 



1 / . 2^ d^ 



1 



d 



1 



2-'-M,2 2^^d^ + 2(' + ^)l + A,2 



With the change of variable |40j 

,2N-1/{2A) 



^'(y,A) = (^(y,A)(l + Ay2)- 
the eigenvalue equation becomes 

(1 + Ay2)/ + (A - 2)y^' + {2£ -1)^ = 
and assuming a power series development 



any ' 



n>0 

the following recursion relation for the coefficients is obtained: 

an 



an+2 = (-1)' 



(n + 2)(n + l) 



[n(An-2 + (2£:-l)]. 
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Therefore the general solution is determined by the values of the coefficients qq £ind ai. In 
particular we can write the solutions yo and yi determined by = l,ai = and gq = 1, ai = 0, 
respectively. The solution yo is an even function while the yi function is odd. The general 
solution will be written as y = aoyo + aiyi. 

The convergence radius of such power series is 



R = — - — because lim 



O'n+2 



|A|. 

This series y = aoyo + aiyi reduces to a polynomial of degree p when one of the two coefficients 
vanish and that means that there exists a positive integer number p such that 2£ — 1 = 2p — Kp^ , 
and then 

The polynomial solutions can be easily found and they have a form quite similar to the corre- 
sponding Hermite polynomials. 



7 Factorization method and shape- invar iance 

The spectrum of the harmonic oscillator can be found using the factorization method we shall 
describe in this section [39] . Actually the Hamiltonian can be factorized and we arrive to a shape- 
invariant Hamiltonian for which the full point-spectrum can be found by algebraic methods. We 
take h = Lo = \. The eigenvalue problem is (up to a factor 1/2) 



where the Hamiltonian H is such that H — 1 
and then the ground state is found from 



i^a, with a = {d/dx + x), = {—d/dx + x), 



aV'o = ( ^ + 2; ) Vo 



0, 



i.e. ipo e ^^Z^. The other eigenstates can be found by applying an iterative way the creation 
operator to the ground state: 

V 2"7i! 

with Hn{x) being the Hermite polynomials. 

A similar procedure can be used for the quantum nonlinear oscillator case. We should look 
for a function W{x), to be called super-potential, in such a way that the operator A and its 
adjoint operator , given by 

A = -L (s/TT^A + w(.) , A- = -L (- VT^.±. + w(.)) . 

are such that Hi = A, i.e. 



A+A = - 
2 



'l + \x'^— + Wix) 
dx 



d 



VT+A^— + W{x) 
dx 
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Therefore, in order to the Hamiltonian so found be that of the deformed quadratic energy term 
together with a potential Vi, the super-potential function W must satisfy 

Vl+Jx^W -W'^ + 2Vi = 0. 
We can define a new quantum Hamiltonian operator 



d 



y^l + Xx^— + W(x) 
ax 



d 



ax 



which is called the partner Hamiltonian. The new potential V2 is given in terms of W by 
V2 = \ [^l + Xx'^W + gW"^^ . 
The important fact is that 

AHi = H2A, A+H2 = HiA+. 

When ^il^-) = El"^), then, #2^1^) = = EAl'i'). If ylj^*) / 0, A]"^) is an eigenvector 

of H2 corresponding to the same eigenvalue E, and similarly, if |$) is an eigenvector of H2 with 
eigenvalue E and such that ^+1$) 7^ 0, then ^"'"l^) is an eigenvector of Hi corresponding to 
the same eigenvalue E. 

The spectra of Hi and H2 are then almost identical, the only differences are when either \^) 
is an eigenvector of Hi but ^|^') 7^ 0, or |<&) is an eigenvector of H2 for which = 0. 

Some parameters may appear in the expression of Vi, and the super-potential function W 
will also depend on them. The most important case is when the explicit forms of the potential 
and its partner are quite similar and only differ in the values of the parameters. In this case we 
say that the problem has shape invariance (see e.g. |41j). 

Suppose that a quantum Hamiltonian Hi{a) admits a factorization Hi{a) = A^{a)A{a) in 
such a way that the partner Hamiltonian H2{a) is of the same form as Hi{a) but for a different 
value of the parameter a. More specifically, there exists a function / such that 

#2(0) = Hi{ai) + R{ai), 

where ai = f{a) and R{a) is a constant depending on the parameter a. In this case Genden- 
shtein developed a method for exact computing of all the spectrum of Hi [lTl[l2lll3]. First, the 
bound state |^o) is found by solving A{a)\^o{a)) = 0, and has a zero energy. Then, |^o(«i)) 
is an eigenstate of H2{a) with Ei = R{ai), because 

^2(a)|^o(ai)) = (^i(ai) + R{ai))\^o{ai)) = i?(ai)|^o(«i)), 

and At(a)|^o(«i)) is the first excited state of ^^1(0), with energy Ei = R{ai), because: 

Hi{a)A\a)\^o{ai)) = A^{a){Hi{ai) + i?(ai))|^'o(«i)) = R{ai)A\a)\^o{ai)). 

Iterating the process we find the sequence of energies for Hi{a) 
k 

Ek = ^R{aj), Eo = 0, 

the corresponding eigenfunctions being 

ao)) = A\ao)A\ai) ■ ■ ■ ^^(Q„_i)|^o(a;, 



where ao = a and aj+i = f{o.j), namely, = f^{ao) = f^{a). 



18 



J.F. Carinena, M.F. Ranada and M. Santander 



Coming back to the nonlinear oscillator case, if /3 G M, we define the linear operator in 



^/2 



f3x 



dx VI + Ax2 J ' 



for which its adjoint operator is 



1 

72 



d 



dx Vl + A2;2 



Then, we find that if H[ = Hi - (l/2)/3, 



I, . 9^ d^ I. d 1 



+ Xx 



H'o = AA+ 



This shows that there is one positive number (3 for which the Hamiltonian H'l of the quantum 
non-linear oscillator admits a factorization with the parameters a and /? related by = /3(/3+A). 
In our preceding case, the parameter being /3, when comparing H'^ with its partner, as 



^1 + ^^'):^ + ^^ 



d 



1 + Ax^ 



1 



(/3-A), 



#((/3-A) = 
then we see that 
H[{13-X)-- 
and therefore 

where / is the function = /3 - A. If i? is defined by R{(3) = /3 + (1/2), then 

This shows that, as the quantum non-linear oscillator is shape invariant, we can develop the me- 
thod sketched before: First, the eigenvector |^o) is determined by the condition A(/3o)|^o) = 0- 
More specifically, we should solve the differential equation 



dx 1 + Ax^ 







and therefore the wave function of the fundamental state must be proportional to 



1 



(l + Ax2)^o' 2A' 
The energies of the first excited states will be 

E[ = R{(3i)=f3-X+'^ 
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and iterating the process we get 



k=l k=l ^ ^ k=l ^ ^ 



and therefore 

,2 

^ - V = n/9 - - 



= n/3 + A 



n/? - —A. 
2 



fc=i 

The energy of the eigenstates of Hi = H[ + (1/2)/? will be given by 

En = nP-—X + -13. 
The method also provides us the corresponding eigenfunctions as 
|^i(/3)) = ^+(/3)|*o(/3i)), 



There is a clear difference between the case A > and the case A < 0. Note that the lowest 
value for E'^ is Eq = 0. Therefore: If A > only values such that 

2 - - A 

are allowed. The eigenvalues are not equally spaced. On the contrary, when A < all natural 
numbers are allowed for n 

8 Quantization of the 2-dimensional nonlinear 
harmonic oscillator 

There exist relatively few examples of quantum Hamiltonians in two dimensions whose spectrum 
can be fully determined by algebraic methods. We shall show that the deformed nonlinear 
oscillator is one of such examples. Furthermore, this can be done in different ways |44j . 

First, one can check that all the measures d/i = p{x,y)dx A dy invariant under the vector 
fields 

d „ r d 



Xi = ^/l + Xr^ — , X2 = VT+Ar^- , 
ox oy 

are proportional to 

dii = , =dx A dy. 

Therefore we shall consider the Hilbert space (M^ , dp.) . 

The form of Killing vectors Xi and X2 generating 'translations' suggests us to take the 
momenta operators 

d ^ r —d 



Px = -Wl + Ar2 — , P = -ihVl + Xr , 
ox oy 
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and then the quantum Hamiltonian is 

H = -—((l + Xr^) — + Xx—] - — ((1 + Xr^) — + Xy—] 
2m \ ' dx^ dx J 2m \ ' dy'^ dy J 



r 



2 



2m \ dy'^ dx'^ dxdy dx dy J 2 1 + Ar^ ' 

which can be written as 



H = Hi + H2- aP, 



with 



^2 = -— 1 + Ar2— - + Ay— +-5 



2m \ dy J 2 1 + Ar^ ' 

^ f od'^ od^ d^ d d\ 1 
J = -t;— [x TTT + y - ^xy^r^ a;- y— + -g 



2m \ dy"^ dx"^ dxdy dx dy J 2 1 + Ar^ 

Note that each term commutes with the sum of the other two and therefore with H. Conse- 
quently, we can consider three different (complete) systems of compatible observables: 

{^1,^2 -AJ2}, {^1-AJ2,#2}, {H1 + H2J}. 

Therefore we should solve one of these spectral problems: 

A) Hi^{Ei,E2j) = Ei^{Ei,E2j), {H2 - XJ^)^{Ei, E2j) = ^2i^(^i, ^2,); 

B) (^1 - Aj2)^(^i„ ^2) = Eij^{E^„E2), H2^{Eij,E2) = E2^{E^j, E2); 

C) Hi^{Ei,E2j) = Ei^{Ei,E2j), {% - XJ^)^{Ei,E2j) = ^2i^(^i, ^2,). 

As in the one-dimensional case, it is convenient to use g = ma^ + Xha and dimensionless 
variables 



h h ma 

X, y = \ — = ^ ^' ^ = ""^ 

ma V ma h 

for which 1 + Xr^ = 1 + Af^ and the Schrodinger equation becomes 

\ ^ dx-^ dx J 2m \ dy^ dy J 

^— 2x x— — -fl A) ''^ -I' - c-I' 

2m \ dy"^ ^ dx"^ ^ dxdy dx ^ dy J 2 1 + Ar^ 

A) As the Hamilton-Jacobi equation separates in coordinates {zx,y) we shall use such coor- 
dinates to write the Schrodinger equation {z is used instead of Zx) 

1 / 1 + Az2 92 d \ ^ 

+ z r^T^ ^ 



2 V 1 + Ay2 9^2 ^y2 

9^ OA C'X, 1 



dy"^ "dzj 2 1 + Ay2 \^ 1 + Az^ 
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and assuming a factorization for ^{z,y) of the form ^{z,y) = Z{z)Y{y) we find the following 
equations: 

-i ((1 + Kz^)Z" + KzZ') + 1(1 + A)^^Z = /.Z, 

- 1 ((1 + WfY" + 2Ay(l + Ky^)Y') + Q(l + A)y2 _ (l + Ay2)e^ y = 

In this way the two-dimensional problem has been decoupled in two 1-dimensional equations. 
The first one is the corresponding one to a 1-dimensional system. The second one however is a 
bit different because includes the contribution of the angular momentum. 

If a new parameter v = e — introduced the second equation becomes 

(1 + A?/2)y" + 2Ayy' - (1 + A - 2A/i)^^^y + 2vY = 0, 
i.e. defining = 1 (1 - 2/i)A, 

(1 + Ky^)Y" + 2Ayy' - gI^^^Y + 2uY = 0, 

and then writing 

y(A,/.) = (/(2/,A)(l + Ay2)-^''/('^\ 
it becomes 

(1 + Ay^)q" + 2(A - G^)yq' + {2u - G^)q = 

which is a deformation of the Hermite equation. 
Assuming the power expansion 

oo 

(Z(2/,A) = ^c4A)y" 

n=0 

we obtain the recursion relation 

^^"+2 = -T—^^r—u (^"(^ - 1) - ^^(2^ + 1) + ■ 

{n + 2)[n + l) 

Therefore the general solution is determined by the values of the coefficients cq and ci. In 
particular we can write the solutions go and qi determined by cq = 1, ci = and cq = 1, ci = 0, 
respectively. The solution q^ is an even function while the qi function is odd. The general 
solution is g = cogo + cigi. 

The radius of convergence of both power series defining q^ and qi is 

1 

R = 



|A| 

The solution reduces to a polynomial of degree n if one of the coefficients is zero and there 
exists an integer number n related with v as follows: 

2u = Gf,{2n + 1) - n{n + 1)A. 

One can study the properties of these polynomial solutions which have a lot of similarities 
with the Hermite polynomials. The associated Sturm-Liouville problem plays a relevant role. 
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In summary, the bound states for this system have energies 
em,n = Hm + I'n = im + n + l) ( 1 - ^(m + n)A 



2 



and wave functions given by 



'I + w 



We finally mention that the properties of these A-dependent Hermite polynomials are dis- 
cussed in [lOl m]. It is proved the orthogonality as well as the existence of a A-dependent 
Rodrigues formula, a generating function and A-dependent recursion relations between polyno- 
mials of different orders. 
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